Emergence of Classical Orbits in Few-Cycle Above-Threshold Ionization 
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The time-dependent Schrodinger equation for atomic hydrogen in few-cycle laser pulses is solved 
numerically. Introducing a positive definite quantum distribution function in energy-position space, 
a straightforward comparison of the numerical ab initio results with classical orbit theory is facili- 
tated. Integration over position space yields directly the photoelectron spectra so that the various 
pathways contributing to a certain energy in the photoelectron spectra can be established in an 
unprecedented direct and transparent way. 

PACS numbers: 32.80.Rm, 42.50.Hz, 34.50.Rk, 02.60.Cb 
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Tunneling ionization in strong laser fields is a prime 
example for nonperturbatively driven quantum systems. 
The complex structure in the photoelectron spectra can 
be interpreted in terms of interfering quantum orbits in 
the spirit of Feynman's path integral approach 0, 0- 
The recently achieved generation of phase-stabilized few- 
cycle laser pulses [jj offers the opportunity to control the 
continuum quantum dynamics of the released electrons 
or, vice verse, to use the electron spectra for measur- 
ing the electromagnetic field of the laser with sub-cycle 
time-resolution 

The hierarchy of theoretical approaches to few-cycle 
above-threshold ionization (FCATI) ranges from "sim- 
ple man's theory" (SMT), considering only the classical 
orbits of the released electrons, via more quantitative, 
semi- analytical theories such as the "strong field approx- 
imation" (SFA) including rescattering of the electron at 
its parent atom and quantum orbit theory (QOT) (see, 
e.g., 0] for a review), to the exact, numerical ab ini- 
tio solution of the time-dependent Schrodinger equation 
(TDSE). The latter yields all observables that can pos- 
sibly be measured in an experiment. However, a real, 
intuitive understanding of the underlying physical mech- 
anisms can only be obtained with the help of simple ap- 
proaches such as SMT. 

This work aims at providing a connection between the 
exact, numerical result and simple physical pictures. To 
that end a quantum distribution function is introduced in 
whose positive definite probability density wave packets 
oriented along classical trajectories emerge as the tun- 
neling ionization regime is approached. The method is 
applied to FCATI where the dominating classical orbits 
are extremely sensitive to the so-called carrier-envelope 
phase (CEP). 

Let \ty(t)) be the numerically determined exact solu- 
tion to the TDSE i9t|*(i)> = H(t)\^{t)) describing the 
outermost electron interacting with a laser field, i.e., 



H{t)= l -[p + A{t)f + V{r) 



(which is well justified for the laser parameters to be 
discussed in this work), and V(r) is the atomic potential. 

The electric field of the laser is given by E[t) = 
—dtA(t). We assume that it is linearly polarized and 
has the form 



E(t) = E(t)e z cos(ujt- 



(2) 



with E(t) being the pulse envelope covering N laser cy- 
cles of period T = 2tt/w, E(t) = E sin 2 [wi/ \2N)] for 
< t < T p = NT and zero otherwise. The CEP (or "ab- 
solute") phase in J2J is irrelevant for long pulses. Since 
recently, however, few-cycle laser pulses with stabilized 
CEP <fi can be generated Under such conditions the 

phase <j) strongly influences the dynamics of the released 
electrons 0- 

The method introduced in this work is based on the 
projection technique proposed in Q. The component 
|$ 7 (£)) of the final wave function |* / ) = |*(T p )) that 
contributes to energies within the bin of width 27 cen- 
tered at £ is calculated as 



\$ 1 (£)) = W 1 {£)\* f ) 



(3) 



where W 7 (£) = j 2 " /[(H - £) 2 " + 7 2 "], and H is the 
Hamiltonian without laser field. With increasing or- 
der n — 1,2,3,... the energy- window W 7 (£) becomes 
more and more rectangular. Numerically, the energy 
component |<& 7 (£)) is calculated by solving the equation 
W 7 T 1 (i?)| < i ) 7 (iS')) = \^f), making use of the factorization 

(H - £) 2 " + 7 2 " = UT=i[H -£ + exp(ii/ Bifc )7][#o - 
£ — exp(it , „ ) fe)7]. The phases v n ,k can be easily calculated 
up to the desired order n. The lowest order values read 
= 7r/2, f2,i — 37r/4, V2.2 — 7r/4. 
The probability P 7 (£) to find the electron in a fi- 
nal state within an energy bin of width 27, order n, 
and centered around £ is P 1 {£) = f \W 2 (£)\^ f ) = 
($ 7 (£)|* 7 (£)) = Z m \(Vf\m)\ 2 /( 7 ,n) where \m) and 
£ m are energy eigenvectors and eigenenergies of Hq, re- 
is 



(1) spectively, and /( 7 , n) = {j 2 " /{(£ m - £) 2 " + 7 2 "]}' 



(atomic units are used throughout). Here, A(t) is the 
vector potential of the laser field in dipole approximation 



the shape-function due to the finite energy-window. 

In our TDSE-solver the wave function is expanded 
in spherical harmonics Y^ m (r2). In the present work 
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we restrict ourselves to linearly polarized laser light 
in dipole approximation so that only to = too with 
Too the magnetic quantum number of the initial state 
contributes, i.e., */(r) = J2e Re( r ) Y e m ° (ty/r- Equa- 
tion © thus yields |$ 7 (£)) in the form $ 7 (£,r) = 

with 

«' 

(4) 

the energy-resolved probability density. P 7 (£,f2,r) may 
be viewed as a quantum distribution function in energy 
and position space. Integration over position space yields 
indeed the electron energy spectrum. However, integra- 
tion over energy yields only approximately the spatial 
probability density. For our purposes, the great advan- 
tage of Q in comparison to, e.g., the Wigner distribu- 
tion function, is (i) that P 7 (£,f2,r) is positive definite 
and (ii) that the energy-width 7 serves as an additional 
parameter, which turns out to play a crucial role for the 
emergence of classical orbits, as will be demonstrated in 
the following. 

Figure H shows the distribution P 7 (£, f2,r) with 7 = 
0.04 for the electron of H(ls) after a 4-cycle pulse with 
CEP = 0, lu = 0.056 (800 nm), and E = 0.0834 
(corresponding to 2.4 x 10 14 Wcm -2 ) in polarization di- 
rection e z (a) and — e z (b). Integration over the posi- 
tion coordinate yields the photoelectron spectra shown 
at the right-hand-side. Note that bound state energies 
£ < are properly treated; the ground state popula- 
tion at £ = —0.5 is clearly visible. Despite the short 
pulse duration, the spectra still show the usual pat- 
tern: after a first plateau up to « 2U P = 1.1 (with 
Up = E 2 /(4:UJ 2 ) the ponderomotive energy), the photo- 
electron yield drops down to the second plateau, reach- 
ing up to w 10Up = 5.5. While for the electrons leav- 
ing in e z direction (a), the latter plateau is well devel- 
oped, the yield continues to decrease more rapidly for 
the electrons in the opposite direction (b). In both di- 
rections a clear interference structure is visible for ener- 
gies in the second plateau. These interferences translate 
to the P 7 -distribution where they appear as horizontal, 
white stripes. In general, the energy-resolved probability 
density shown in Fig.^ is rather delocalized in position 
space (note that more than 1000 a. u. in ^-direction are 
covered). This is expected since the photoelectrons con- 
tributing to a certain energy may be well described by 
plane waves once the laser pulse is over and they are suffi- 
ciently far away from the origin. These plane waves may 
interfere destructively, giving rise to the complex pattern 
visible in the photoelectron spectra. The rather sharp 
interference patterns in the spectra of Fig.^ are thus al- 
ways accompanied by a spatially delocalized density P 7 . 
SMT, on the other hand, deals with classical trajectories 
of point-like electrons. 
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FIG. 1: Logarithmically scaled contour plots of P 7 - 
distributions (7 = 0.04, n = 3) after a 4-cycle pulse with 
<j> = 0, w = 0.056, E = 0.0834 in (a) e z -direction, i.e., 
Q, = (1?, ip) = (0, 0) (ip is arbitrary for linear polarization 
along e z ), and in (b) —e z -direction fl = (71", 0). Directional 
photoelectron spectra at the right-hand-side are obtained by 
integrating P 1 over r with i? fixed. 

Let us briefly review SMT: an electron is "born" at 
time to with zero initial velocity and from then on moves 
in the laser field as if there were no Coulomb potential. 
In the case of a laser field polarized in z-direction, the 
velocity at times t > t Q is vi°\t) = — f*dt' E(t') = 
A(t) — A(to) so that the final photoelectron energy and 
the position are given by 

^l = l[MT P )-A(t )} 2 (5) 

and 

Z W(t) = z(t)-z(to)-A(to)(t-t ) (6) 

with z(t) = Jgdi' A(t'), respectively. 

The upper index (0) in and © indicates the "di- 
rect" electrons that do not interact anymore with the ion 
after their emission. If, instead, the electron is driven 
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back to the origin, it may rescatter at time t r if the 
condition |^^(t r )| < d r is fulfilled. The smaller the 
distance d r is allowed to be, the less rescattering solu- 
tions exist in SMT. Upon rescattering electrons may as- 
sume a lateral velocity. Let X denote the scattering angle 
with respect to the velocity vector of the incoming elec- 
tron. The velocity after the (elastic) scattering event 
then reads (using cylindrical coordinates z,p) v^(t) — 
[vf\t r )cosx + A(t) - A(tr)]e z + \vi"\t r )\ sm X e p with 
v^z\t r ) = A(t r ) — A(to). For the expected electron posi- 
tions 

zW(t) = z(t)-z(t )-A(t )(U-t ) (7) 

+[A(t r )(co8x-l)-A(t )cosx\{t-t r ), 
p<U(t) = [A(t r )-A(to)}{t-t r )sm X (8) 

follow, and the final energy is given by = 
v^(T p ) 2 /2. In the extreme case of 180°-deflection x — 77 
one has £^.(x=*) = [A(T p ) + A(t Q ) - 2A{t r )f /2. 

SMT is able to predict all the possible positions 
z (0) (T p ) and z^(T p ), p (1) (T p ) of an electron with final 
energy £ = £^ = £^. The higher is the ionization 
probability at time f (and the rescattering probability 
at time t r ) the higher should the probability to actually 
measure an electron at the positions predicted by SMT 
be. In QOT, this is taken into account on the level of 
an extended SFA transition amplitude with rescattering 
included @. 

Let us now compare the numerically determined P 7 - 
distribution with the SMT predictions. In Fig.[2]the final 
positions z^ (T p ) and z^ (T p ) are presented as a function 
of the emission time to for the laser pulse with parame- 
ters as given in the caption of Fig.Q] The final positions 
z(°)(T p ) (direct electrons) are plotted in bluish color, the 
final positions z^'(T p ) of the rescattered electrons with 
X = 7r arc plotted reddish. The lighter the color is, the 
higher was the absolute value of the electric field am- 
plitude at the time of emission (i.e., the higher was the 
ionization probability). 

Figure [3] shows the P 7 -distributions corresponding to 
Fig-E now calculated with a ten times wider energy- 
window (7 = 0.4). The SMT predictions of Fig. are 
included using the same color coding. The larger energy 
width allows for rather localized electron wave packets in 
position space because of the uncertainty principle with 
respect to P and r. The wave packets are clearly ori- 
ented along the classically expected positions for suffi- 
ciently high energies. The dotted spectra on the right- 
hand-side were calculated with 7 = 0.4 as well, clearly 
showing that the emergence of the classical orbits in the 
energy-resolved probability density P 7 is accompanied by 
a loss of energy resolution. Compared to the solid spec- 
tra (calculated with 7 = 0.04), the interference pattern 
in the rescattering plateau is completely washed out in 
(a) and barely visible in (b). Instead, a constructive in- 
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FIG. 2: (Color online) Final positions z m (T p ) and z w (T p ) 
vs emission time to for direct electrons (d) and rescattered 
electrons with \ = it (rl-r5). The shape of the laser field is in- 
dicated (dashed). The rescattering condition was Iz*- ' (i r )| < 
d r = 1. Only emission times to where |i?(fo)| > E/10 were 
taken into account in the calculation of z"'(Tp). The lighter 
the color is, the higher is the absolute value of the field ampli- 
tude at the time of emission. Same color-coding and labelling 
are used in Fig.[!|] 



terference between the rescattering branches r2 and r4 is 
observed in the P 7 -distribution of Fig.|U>. 

The result for the electrons in — e^-direction in Fig.^ 
is particularly interesting because of the two competing 
rescattering solutions reaching up to high energies. As 
expected from SMT, the dominating solution first fol- 
lows the more probable direct pathway (the lighter blue 
branch up to \z\ « 400). At energies beyond the direct 
cut-off £ w 1.1, the quantum mechanical result continues 
to follow this branch despite the fact that the rescat- 
tering branches r3 and r4 are more likely. This is be- 
cause, contrary to the sharp cut-offs in SMT, the transi- 
tion from the direct plateau in the spectrum down to the 
six orders of magnitude less probable rescattering regime 
is smooth in the quantum mechanical result. Once the 
level of the rescattering plateau is reached, the dominat- 
ing branch is given by the SMT solution r4. The energy 
cut-off of r4 is again quantum mechanically extrapolated 
up to £ ps 4.5 before the less likely but to higher energy 
extending rescattering solution r2 overtakes. It is this 
branch r2 that defines the highest possible classical en- 
ergy £ « 5 and thus defines the ultimate classical cut-off 
for the electrons in — e 2 -direction. The preceding analysis 
nicely demonstrates the usefulness of the P 7 -distribution: 
through the comparison with the expected electron po- 
sitions in SMT, the contribution of various pathways to 
each energy in the exact, numerically determined pho- 
toelectron spectrum can be understood in simple terms 
and great detail. 

So far our study was restricted to electrons moving 
along the laser polarization direction. Figure 01 shows 
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FIG. 3: (Color online) Emergence of electron wave packets 
in the P 7 -distributions with 7 = 0.4, n = 3 for the same laser 
pulse as in Fig.0 Final positions as predicted by SMT (cf. 
Fig.|2J are included using the same labels and color-coding. 
Dotted spectra at the right-hand-side were also calculated 
with 7 = 0.4. 

P 7 -distributions (7 = 0.04) at time t = T p in the pz- 
plane for fixed energy £ = 5.0 and different CEPs tj> — 
0, tt/4, tt/2, and 3vr/4. The expected SMT positions 
according J7J) and (JSJ are included in the plots. Note 
that for the first three phases <fi no SMT solutions exist 
for electrons traveling in — e z -direction since £ — 5.0 lies 
beyond the corresponding classical cut-off in these cases. 
The classical solutions (where they exist) agree well with 
the quantum mechanical results. At high energies, the 
electron emission into narrow cones oriented along the 
laser polarization direction is predicted by SMT and well 
confirmed by the numerical result. Figure 0] illustrates 
the strong phase-dependence of few-cycle ATI. For <f> — 
7r/4 almost no electrons of the selected energy are emitted 
in — e z direction. For tf> = 7r/4 and 7r/2 a destructive 
interference for the emission angle D = 0.245 (dotted line 
in (b) and (c)) is observed, which is absent for the other 
two phases where maximum emission of electrons with 
£ = 5.0 is not in e z -direction but slightly off-axis. 

In summary, a positive definite quantum distribution 
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FIG. 4: (Color online) P 7 -distributions (7 = 0.04, n = 3) at 
t — T p in the pz-plane for £ — 5.0 and (a) = 0, (b) 7r/4, 
(c) 7r/2, (d) 37t/4. SMT solutions are included. Dotted line 
indicates destructive interference for 1? = 0.245. 



function P 7 with variable energy resolution 7 for the iden- 
tification of relevant quantum orbits has been introduced. 
Upon decreasing the energy resolution, the emergence of 
localized wave packets is observed if the phenomenon un- 
der investigation is accessible to (quasi) classical theory. 
The proposed method has been applied to the numeri- 
cally determined, exact electron wave function for H(ls) 
in few-cycle laser pulses. Pronounced interference pat- 
terns at high energies have been observed in the photo- 
electron spectra. Using the P 7 -distribution, the complex 
structure of the electron spectra has been analyzed in 
terms of interfering, delocalized electrons at high energy 
resolution and localized electron wave packets at lower 
energy resolution. In that way the connection between 
the classical orbits of so-called "simple man's theory" and 
quantum orbit theory has been established, starting from 
the exact, numerically determined wave function. 

The author thanks Dejan Milosevic and Wilhclm 
Becker for illuminating discussions. This work was sup- 
ported by the Deutsche Forschungsgemeinschaft through 
a Heisenberg fellowship. 
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